Shot noise in frustrated single-electron arrays 
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We have carried out numerical simulations of shot noise in 
2D arrays of single-electron islands with random background 
charges. The results show that in contrast with the ID arrays, 
at low currents the current noise is strongly colored, and its 
spectral density levels off at very low frequencies. The Fano 
factor may be much larger than unity, due to the remnants 
of single-electron/hole avalanches. However, even very small 
thermal fluctuations reduce the Fano factor below 1 for almost 
any bias. 

PACS numbers: 73.23.Hk, 73.40.Rw, 85.35.Gv 

Arrays of small conducting islands separated by tun- 
nel junctions are among the simplest systems exhibiting 
the Coulomb blockade and the related effects of corre- 
lated single-electron tunneling - for reviews see, e.g., Refs. 
[1-3]. Besides the successful use of multi-junction arrays 
for absolute thermometry (see, e.g., Ref. [4]), two their 
important applications in future digital nanoelectronics 
have been suggested: 

(i) The arrays may replace the double junctions in 
single-electron transistors, in order to raise the single- 
electron addition energy and hence increase the opera- 
tion temperature range [5]. A quantitative analysis of 
this opportunity [6] has shown that for this purpose, ar- 
rays with zero background charges of the islands may be 
especially effective. 

(ii) The arrays may serve as circuit elements with " sub- 
electron" (quasi-continuous) charge transfer [7], capable 
of leaking to ground the random background charges of 
single-electron islands induced by charged impurities - 
see, e.g., Ref. [8]. Earlier we showed [9] that both condi- 
tions of sub-electron transfer, that had been formulated 
in the introduction to Ref. [7] , are satisfied if the back- 
ground charges of the array islands are arranged in a 
special way. 

Both applications are, however, strongly affected by 
the background charge randomness. In our recent work 
[10] we have shown that in the realistic case of "full 
charge frustration" , when random charges are spread uni- 
formly within the available range [— e/2, +e/2] [11], the 
former application becomes impossible, because the ar- 
ray loses any sensitivity to the gate voltage. Concerning 
the second application, the results of the analysis [10] 
were somewhat more encouraging: they showed that the 
background charge randomness suppresses the Coulomb 
blockade considerably, albeit not fully. Moreover, in 
large 2D arrays the statistical distribution of the block- 
age threshold voltage Vt becomes relatively narrow, so 
that dynamic properties of these arrays are virtually re- 
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FIG. 1. Spectral density of current noise calculated for sev- 
eral 2D arrays with random background charge distribution, 
as a function of observation frequency u), for two values of 
dc current /. For / = O.le/RC only two extreme cases are 
shown. The inset shows the corresponding results for ID ar- 
rays 



producible sample-to-sample. In the light of this result, 
it was important to check whether the second condition 
of sub-electron transfer [7] , the shot noise suppression, is 
also satisfied in such arrays. 

Unfortunately, despite the recent successful experi- 
mental implementation of ordered 2D arrays, whose 
transport properties are well described by the full- 
frustration model, by several groups (see, e.g., Refs. 12- 
16), we are not aware of any measurements of their noise 
properties. In this situation, the objective of this work 
was a numerical simulation of shot noise in 2D arrays 
with random background charge distribution, within the 
framework of the "orthodox" theory of single-electron 
tunneling [1]. For comparison, results for ID arrays are 
also shown. (Some preliminary results for the ID case 
were published earlier [7].) For the sake of simplicity we 
have studied the most interesting case of small island self- 
capacitance: Co <C C/N 2 , where C is the capacitance be- 
tween neighboring islands, and N is the array size. The 
2D arrays were of an almost square shape, with N — 1 
islands along the current direction and N islands across. 
Our algorithms were based on the well-known Monte 
Carlo code MOSES 1.2 that is available on the Web at 
http://hana.physics.sunysb.edu/set/software/index.html. 
Despite the code optimization, getting reasonable accu- 
racy of noise calculations has required substantial super- 
computer resources. For example, the calculation of just 
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FIG. 2. The Fano factor F = Si(0)/2el as a function of dc 
current per row for charge-frustrated arrays of two different 
sizes, and for the 20 x 20 array with zero background charge 
[20]. The inset shows the corresponding plots for frustrated 
ID arrays. 
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FIG. 3. The average "crossover" current I cr (the dc cur- 
rent corresponding to F = 1/N 1 ^ 2 ) as a function the array 
length for ID and 2D arrays, both with and without charge 
frustration. Error bars correspond to the r.m.s. width of the 
statistical distribution of I cr . 



one group of curves shown in Fig. 2 took approximately 
1,500 node-hours of modern CPUs. 

Figure 1 shows typical frequency dependences of the 
spectral density Si(iv) of current noise in 2D arrays, for 
two values of dc current / (determined by dc voltage 
V > V t applied to the array) at zero temperature. At 
large currents and/or high frequencies the results are 
virtually independent on the system dimensionality; in 
fact, general arguments [7,17] show that in all cases S/(w) 
tends to 2eI/N. However, at finite currents the results 
for 2D and ID cases are quite different. In contrast to the 
ID case, 2D arrays typically do not exhibit the spectral 
density increase at uj/2tt = I/e, corresponding to the 
so-called SET oscillations [1]. More importantly, as dc 
current is reduced and/or the array size is increased, the 
frequency at which S7(w) levels off decreases rather dra- 
matically [18], while the saturation level 57(0) increases 
well beyond the Schottky formula value 2el. (In ID ar- 
rays, 5*7(0) is always below 2el, so that the Fano param- 
eter F = S[(Q)/2eI is always below unity [7].) Figure 2 
shows the Fano factor as a function of dc current for sev- 
eral arrays with random background charge, and for one 
20 x 20 array with zero background charge. At J — ► 0, 
the Fano factor saturates for any particular system, but 
at a very high level, with exponentially broad statistical 
distribution. 

Note that the dramatic increase of noise in 2D arrays 
is limited to low currents. This is why it does not affect 
too much such measure of the shot noise suppression as 
the "crossover" value I cr of dc current (Fig. 3), at which 
the Fano factor drops to 1/iV 1 / 2 , i.e. to the geometrical 
mean between the unity and its asymptotic, high-current 
value [7]. 

The difference between results for 2D and ID arrays 
may be attributed to the different topology of charge car- 



rier motion. In the ID case, current is carried along the 
array by single stream of single electrons (holes). With- 
out disorder, the array provides a single bottleneck for 
this motion. At / — > 0, overcoming this bottleneck be- 
comes a dominating source of current fluctuations, so 
that the current fluctuations are essentially the same as 
in a single tunnel junction, with a broadband shot noise 
and F — ► 1 [7]. Disorder in a long array induces mul- 
tiple quasi-bottleneck states with very close thresholds 
V} 1 \ The actual threshold voltage V t is determined by 
the highest of these thresholds, so that at V — > Vt (I — ► 0) 
the Fano factor still tends to unity. 

On the contrary, broad 2D arrays may feature mul- 
tiple streams of electrons and holes moving in opposite 
directions, that do not necessarily annihilate even in the 
absence of disorder, giving rise to long "avalanches" pro- 
ducing in particular very large noise (F 3> 1) [19]. The 
disorder suppresses this effect (Fig. 2), but only partly. 

It is interesting (and important for applications) that 
the Fano factor at low currents is substantially sup- 
pressed by even very small thermal fluctuations. For 
example, Fig. 4 shows that temperature as low as 
ksT = 10~ 2 e 2 /C reduces the Fano factor at low cur- 
rents by several orders of magnitude, bringing it below 
unity, besides a very small region where eV < 2ksT. (In 
that region, the current fluctuations are essentially ther- 
mal, rather than due to the shot noise.) This behavior is 
very much different from that of ID arrays (inset) where 
small increase of temperature have negligible effect on 
shot noise. 

Actually, Fig. 4 may be somewhat misleading, because 
at vanishing temperature the limit / — > corresponds to 
the voltage approach to a finite value V t , while at T > 
it is achieved at V — > 0, so we are speaking about very 
different dc bias regions. In fact the decrease of F is 
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FIG. 4. Fano factor as a function of dc current for sev- 
eral 20 x 20 arrays, at zero and finite temperatures. The inset 
shows the corresponding dependences for a typical charge dis- 
tribution of a ID array. 



not a result of the decrease of current fluctuations as 
such: actually at fixed applied voltage the spectral den- 
sity 5*7(0) always grows with temperature, but dc current 
grows even faster, resulting in the reduction of F. Nev- 
ertheless, Fig. 4 relays the important information, since 
it is the Fano factor, not 5/(0) per se, that is important 
for applications [7]. 

In fact, the reduction of F by thermal fluctuations has 
been already noted in our work on avalanches in uni- 
form arrays (see Fig. 5 of Ref. [19]) and interpreted 
as avalanche suppression due to larger side diffusion of 
electrons and holes, increasing chances for their recombi- 
nation. Apparently, disorder enhances this process and 
suppresses avalanches even further. 

To summarize, we have found that large 2D single- 
electron arrays biased close to the Coulomb blockade 
threshold at T — > (and hence carrying very low cur- 
rent) exhibit strongly colored noise, following approxi- 
mately the l/f a dependence until very low frequencies 
where it levels off yielding very high values of the Fano 
factor. This effect, that may be interpreted as a result 
of electron/hole avalanches in 2D arrays [19], would be 
very unfavorable for the application of such arrays as cir- 
cuit components with quasi-continuous ("sub-electron") 
charge transfer. However, at small but Unite tempera- 
tures the Fano factor is suppressed below unity for vir- 
tually all bias voltages. 
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